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ABSTRACT 

Multiscale topological complexity of the solar magnetic field is among the primary factors controlling energy 
release in the corona, including associated processes in the photospheric and chromospheric boundaries. We present 
a new approach for analyzing multiscale behavior of the photospheric magnetic flux underlying these dynamics as 
depicted by a sequence of high-resolution solar magnetograms. The approach involves two basic processing steps: 
(1) identification of timing and location of magnetic flux origin and demise events (as defined by DeForest et al.) 
by tracking spatiotemporal evolution of unipolar and bipolar photospheric regions, and (2) analysis of collective 
behavior of the detected magnetic events using a generalized version of the Grassberger-Procaccia correlation 
integral algorithm. The scale-free nature of the developed algorithms makes it possible to characterize the dynamics 
of the photospheric network across a wide range of distances and relaxation times. Three types of photospheric 
conditions are considered to test the method: a quiet photosphere, a solar active region (NOAA 10365) in a quiescent 
non-flaring state, and the same active region during a period of M-class flares. The results obtained show (1) the 
presence of a topologically complex asymmetrically fragmented magnetic network in the quiet photosphere driven 
by meso- and supergranulation, (2) the formation of non-potential magnetic structures with complex polarity 
separation lines inside the active region, and (3) statistical signatures of canceling bipolar magnetic structures 
coinciding with flaring activity in the active region. Each of these effects can represent an unstable magnetic 
configuration acting as an energy source for coronal dissipation and heating. 
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turbulence 
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1. INTRODUCTION 

The solar magnetic field is often seen as the underlying cause 
of solar activity and the primary agent coupling the Sun to the 
heliosphere. High-resolution observations show that it possesses 
an incredible level of complexity across all relevant spatial 
and temporal scales (e.g., Schrijver et al. 1992; Balke et al. 
1993; Lawrence et al. 1993; Cadavid et al. 1994; Abramenko 
et al. 2002; Abramenko & Yurchyshyn 2010a; Georgoulis 2005; 
McAteer et al. 2010; Parnell et al. 2009). Understanding the 
origin and the physical significance of this complexity is among 
the most fundamental unresolved problems of contemporary 
solar science. Its further analysis requires accurate quantification 
of cross-scale plasma processes involved in the formation and 
evolution of the observed multiscale magnetic patterns. 

Solar magnetic complexity is deeply rooted in the life cycle 
of coronal perturbations. The coronal magnetic field is known 
to accumulate progressively large amounts of magnetic energy 
and helicity before producing explosive phenomena such as 
flares and coronal mass ejections (see, e.g., Solanki et al. 2006). 
The emergence of complex domains of magnetic connectivity 
accompanying this process is instrumental as it makes the 
magnetic free energy available for dissipation (Antiochos 1998; 
Schrijver et al. 2005; Edmondson et al. 2010). The initial buildup 
of the magnetic complexity can be driven by the injection of 
new magnetic structures (e.g., Falconer et al. 2008; Abramenko 
& Yurchyshyn 2010a; McAteer et al. 2010; Conlon et al. 
2010), or by the fragmentation of the existing flux subject to 
turbulent footpoint shifting both in solar active regions (see, 
e.g., Abramenko & Longcope 2005) and in the quiet Sun (Lamb 
et al. 2008; Chaouche et al. 201 1 and references therein). 


The resulting magnetic patterns featuring mixed-polarities 
regions with strongly sheared fields are prone to kinetic insta- 
bilities which can unfreeze the magnetic flux and start magnetic 
reconnection. The reconnection brings the system to a lower 
energy state and reduces the magnetic complexity locally but 
usually not globally (Aschwanden 2006). Cooperative interac- 
tions between multiple reconnecting plasma tubes may lead to 
the second generation of multiscale phenomena at the coro- 
nal level, including self-organized criticality (Lu & Hamilton 
1991), storms of nanoflares in individual multistranded coronal 
loops (Klimchuk 2006; Viall & Klimchuk 2011; Morales & 
Charbonneau 2008), and cascades of instabilities involving mul- 
tiple reconnecting loops (Hughes et al. 2003; Torok et al. 201 1). 
The distortions in the topology of open field regions caused by 
some of these effects can be a major factor controlling both 
slowly varying and transient solar wind (Antiochos et al. 2007). 
In the quiet Sun, the wide ranges of spatial and temporal scales 
manifested by the brightness of the soft X-ray and EUV corona 
are likely to reflect the multiscale geometry and dynamics of 
the underlying photospheric network (Falconer et al. 1998). The 
importance of solar magnetic complexity in triggering coronal 
dissipation events and the possibility of their self-organization 
into multiscale structures is confirmed by many theoretical mod- 
els (see, for instance, Charbonneau et al. 2001; Aschwanden & 
Parnell 2002; Aschwanden 2011). 

Gaining a deeper insight into the nature of solar mag- 
netic complexity and its impact on coronal instabilities 
requires innovative data analysis tools focused on cause 
and effect links between two or more interconnected so- 
lar processes. There is a serious need to explore multi- 
scale cross-correlations between these processes, in addition 
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to the autocorrelation properties studied before. Furthermore, 
the causality of solar activity is an inherently spatiotemporal 
phenomenon which unfolds in spatial and temporal domains of 
analysis. Until now, the spatiotemporal aspect of solar magnetic 
complexity has not been given proper attention in observational 
studies. Thus, for example, the algorithms of solar magnetic 
tracking commonly used for analyzing the emergence, cancel- 
lation (Martin 1984), fragmentation and coalescence dynamics 
of the photospheric network are essentially monoscale. The out- 
puts of these algorithms and the interpretation of the obtained 
results depend, explicitly or implicitly, on ad hoc choices of 
spatial and temporal correlation scales (see DeForest et al. 2007 
for a detailed review). To our knowledge, none of the methods 
proposed so far can simultaneously address short- and long- 
range properties of the solar magnetic network within the same 
data analysis framework while treating both spatial and tem- 
poral scales as free parameters. The available post-processing 
techniques, such as Fourier spectrograms, spatial correlation 
functions, or multifractal statistics, share the same methodolog- 
ical “blind spot” — they cannot adequately describe transient 
cross-scale coupling associated with the time-dependent nature 
of studied multiscale patterns. 

In this paper, we propose a new data analysis paradigm 
focused on multiscale spatiotemporal properties of the pho- 
tospheric flux dynamics. As opposed to previously proposed 
methods for identifying correlated magnetic events, the devel- 
oped algorithms contain no a priori assumptions regarding the 
intrinsic spatial and temporal interaction scales. By accepting 
this strategy, we take into account the fact that magnetic con- 
centration regions, acting as nodes of the photospheric network, 
can communicate over arbitrarily long distances beyond the lo- 
cal neighborhood. We also assume that these communications 
can encompass a wide range of response times. 

Our technique involves two main processing steps. First, we 
determine timings and locations of origin and demise magnetic 
events (DeForest et al. 2007) based on direct spatiotemporal 
tracking of unipolar and bipolar photospheric regions above a 
specified magnetic flux threshold. Second, we analyze pairwise 
correlations between the detected magnetic events using a gen- 
eralized Grassberger-Procaccia correlation integral technique. 
By design of the algorithms, any detected characteristic scale is 
an attribute of the studied data and not that of the method. 

The proposed data analysis framework is successfully tested 
on a sequence of high cadence Solar and Heliospheric Obser- 
vatory Michelson Doppler Imager ( SOHO MDI) magnetograms 
representing a solar active region (NOAA 10365) and a quiet 
Sun. The application of spatiotemporal methodology extends 
earlier statistical analyses of multiscale geometry of the photo- 
spheric field (e.g., Lawrence et al. 1993; McAteer et al. 2010), 
and allows us to evaluate scaling characteristics of the dynam- 
ics of the photospheric flux not reported before. The results 
obtained reveal a significant difference in the network dynam- 
ics underlying flaring and non-flaring regimes of NOAA 10365, 
and deliver a detailed picture of spatiotemporal clustering in the 
adjacent quiet photosphere. 

The paper has the following structure. Section 2 provides the 
description of the studied data, solar conditions, and applied 
statistical tools. Section 3 presents the main results of our 
analysis and their qualitative interpretation in terms of key 
multiscale plasma processes operating in the solar photosphere 
and corona. Section 4 summarizes our findings and highlights 
some remaining problems. 


2. DATA AND METHODS 
2.1. Observational Data and Solar Conditions 

We analyzed a sequence of 816 high-resolution SOHO MDI 
magnetograms (Scherrer et al. 1995) collected between 14:00 
05/26/2003 and 04:35 05/27/2003, with an average cadence 
time of 64 seconds. The images were calibrated using standard 
SolarSoft procedures and rebinned down to the 0.9 x 0.9 Mm 
resolution to reduce the pixel noise, with a final image size 
of 512 x 512 pixels. The studied field of view contains two 
adjacent spatial domains representing different levels of solar 
activity (Figure 1). 

The quiet region to the north of the solar equator (labeled 
with letter Q in Figure 1) exhibits a characteristic pattern of 
strongly fragmented magnetic flux loosely organized into an 
irregular cellular structure typical of a quiet Sun. As we dis- 
cuss below, this structure is associated with photospheric meso- 
and supergranulation, and it demonstrates a significantly geo- 
metric asymmetry of positive and negative magnetic polarities 
precluding direct identification of canceling bipolar pairs. The 
region operates in a statistically steady state as evidenced by the 
similarity of scaling properties of origin and demise dynamics 
across the studied range of scales. 

The southern portion of the studied field of view (labeled with 
A in the figure) is dominated by the NOAA 10365 active region 
which first appeared on the visible solar disk on 2003 May 19 
at heliographic latitude 5°. A few days later, a new bipolar 
structure formed inside the active region. By merging with the 
decayed preexisting flux, the new dipole developed a complex 
polarity separation line (Verbeeck et al. 2011). Starting from 
25 May, the total magnetic flux of NOAA 10365 was increasing 
steadily before reaching a plateau on 2003 May 29. The region 
gradually rotated counterclockwise with the leading negative 
flux polarity (Tian et al. 2011). The continuous injection of 
a new flux in the northern portion of NOAA 10356 created 
a S — spot configuration with a positive region progressively 
surrounded by the emergent negative flux (Chandra et al. 2009). 
The magnetic field pattern of the emerging bipole featured two 
large elongated areas of opposite polarity, suggesting a twisted 
flux tube of positive helicity (Tian et al. 2011; Chae et al. 2004). 

The activity of the NOAA 10365 region underwent significant 
changes during the studied time interval. Between 15:00 05/ 
26/2003 and 20:48 05/26/2003, the active region generated 
multiple M-class flares, with the GOES 12 1. 0-8.0 A flux 
remaining above the 10” 10 W m -2 level. During 20:48 05/26/ 
2003-02:42 05/27/2003, the GOES flux stayed consistently 
below this level, indicating no noticeable flare activity. The two 
intervals were analyzed separately in order to compare magnetic 
clustering effects in flaring and non-flaring regimes. 

2.2. Definition of Magnetic Events 

To determine timings and positions of magnetic events in the 
studied set of magnetograms, we invoked a direct spatiotem- 
poral tracking approach (Uritsky et al. 2007, 2010a) in which 
magnetic features living for more than one image sampling in- 
terval were treated as subvolumes of the three-dimensional (3D) 
space-time. The history of each feature is represented by a con- 
nected set of image pixels containing an unsigned magnetic field 
above a specified detection threshold B t h- 

The idea of the threshold-based spatiotemporal tracking has 
been implemented in the context of solar magnetograms by 
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Figure 1 . SOHO MDI magnetogram (14:00 05/26/2003) showing the studied field of view (approx. 445 x 445 Mm, centered at the solar disk center) including the 
subregions with quiet (active) conditions above (below) the dashed horizontal line representing the solar equator. The NOAA 10365 active region is outlined by the 
dotted circular line. 


Berger & Title (1996) and Berger et al. (1998). The object 
tracking method used in their studies was based on measuring 
the centroid locations of the most compact and intense magnetic 
features. Here we invoke a somewhat different approach which 
does not rely on the centroid location as a guiding center for 
feature tracking. Time-evolving magnetic concentrations are 
treated as true 3D objects, which makes it possible to study 
a more general class of magnetic structures whose positions 
cannot be adequately described by a centroid. 

In order to overcome the inherent memory limitations of stan- 
dard cluster analysis algorithms such as, e.g., the Label Region 
function of the Interactive Data Language (IDL), we have de- 
veloped a new numerical scheme enabling fast decomposi- 
tion of a multidimensional data array into a set of connected 
pixels (clusters) while dramatically reducing the amount of 
stored information. 

The first step of our technique consists of building a table of 
so-called “activations” — contiguous time intervals during which 
the condition |B| > B lh is fulfilled in each of the image pixels. 
The second step is to find and label spatially connected clusters 
of the activations. This is done using the “breadth-first search” 
principle which avoids backtracking of search trees representing 
individual clusters. We find that it is important to consider all 
of the nearest neighbors of each spatiotemporal node, including 
the diagonal neighbors, when identifying connected activations. 
Finally, the table of activations is rearranged according to the 
assigned cluster labels, in order to provide faster access to the 
detected structures. The output data array contains the complete 
information on the location, shape, and time evolution of all of 


the detected magnetic concentration regions in the studied 3D 
volume. 

It can be seen that the proposed algorithmic solution reverses 
the usual order of feature tracking steps (feature segmentation 
within in an individual image followed by cross-frame feature 
identification within an image set, see, e.g., DeForest et al. 
2007). Instead, we start with detecting temporal traces of 
features in individual photospheric locations, and then consider 
their spatial adjacency. We found this algorithmic solution to 
be remarkably fast and memory-efficient, with the size of the 
output data array being typically two orders of magnitude 
smaller compared to the original data. The proposed scheme 
is obviously easy to parallelize. A more complete description 
of the developed algorithm and its numerical validation on 
high-Reynolds number magnetohydrodynamic turbulence can 
be found in Uritsky et al. (2010b). 

Using the described tracking technique, we identified spatial 
and temporal coordinates of origin and demise magnetic events 
described by positive and negative line-of-sight magnetic polar- 
ity. The origin events (OEs) represent newly detected magnetic 
concentration regions which just exceeded the threshold. The 
demise events (DEs) occur when a previously found region be- 
comes undetectable, i.e., when its flux content drops below B t \ t . 
This classification is analogous to the subdivision into “origin” 
and “demise” events described by DeForest et al. (2007) al- 
though these authors used a different feature tracking method. It 
should be emphasized that our definition of events is not limited 
to the injection of new magnetic bipoles into the photosphere 
and their subsequent submergence into the convection zone. 


3 





The Astrophysical Journal, 748:60 (12pp), 2012 March 20 


Uritsky & Davila 


In essence, we deal with an apparent dynamic associated with 
crossings of the detection threshold B± by the local magnetic 
field. While the resulting sets of events may contain a fraction 
of truly emerging and submerging magnetic bipoles, they also 
represent magnetic flux coalescence, fragmentation, cancella- 
tion, and other effects leading to spatiotemporal variability of 
the photospheric field. 

The photospheric velocities of positive and negative mag- 
netic spots inside the NOAA 10365 active region were in the 
range 0.07-0.24 km s -1 (Chandra et al. 2009), or less than 
0.015 Mm minute -1 . For the timescales r < 20 minutes stud- 
ied here, this motion would produce a spatial displacement of 
^0.2 Mm, which is below the spatial resolution of SOHO MDI 
data. Based on this estimate, the magnetic flux configuration 
of the active region can be considered quasi-static in the pho- 
tospheric plane, with most of the observed changes owing to 
enhancement and reduction of local photospheric magnetic flux 
rather than its horizontal convective motion. The same conclu- 
sion applies to the quiet region. Typically, the quiet-Sun hori- 
zontal speeds lie in the range 0.45-0.50 km s -1 (Shine et al. 
2000); over a 20 minute interval, this motion is likely to be 
limited to a single MDI pixel and remain unresolved. 

To reduce the level of granulation noise (DeForest et al. 
2007), we focused on the magnetic features which lasted more 
than 5 minutes and encompassed a sufficiently large spatial 
area (>4 Mm 2 ). The valid features were also required to have 
small initial and final linear sizes (<(3 Mm) ensuring accurate 
determination of their locations. The magnetic features which 
did not meet these quality constraints, as well as those truncated 
by the boundaries of the chosen field of view or observing 
interval, were excluded from further analysis. 

2.3. The Generalized Correlation Integral Analysis 

The structure of multiscale data sets such as the ones 
generated by solar magnetograms can be accurately described 
using the family of ordcr-c/ generalized fractal dimensions 
(Grassberger 1983): 

' a E^iV/iogp.- 

n 3 log € 

" 1 3 log E "=1 pi 

q — 1 3 log e 

Here, n(e ) is the number of non-empty hyper-cubes of a 
partitioning with length scale e, and p, is the probability for 
a data point to fall into a cube i of the partitioning; see Schuster 
& Wolfram (2005) for a pedagogical introduction to the concept. 
Among the dimensions defined by Equation (1), the Hausdorff 
fractal dimension ( q — 0), the information dimension (q = 1), 
and the correlation dimension (q = 2) are widely used. We 
base our analysis on the latter because it represents a second- 
order statistical moment which can be naturally converted into 
a cross-correlation measure of bipolar magnetic structures, as 
shown below. 

For a set of point objects embedded in a multidimen- 
sional metric space, the correlation dimension is commonly 
evaluated using the Grassberger-Procaccia correlation integral 
(Grassberger & Procaccia 1983). The correlation integral (Cl) 
is computed by counting the number C of pairs of set elements 
separated by the Euclidean distance < i. The shape of the C(l) 
dependence characterizes clustering effects in the studied set. 
For a random set with statistically independent element loca- 
tions C(l ) oc l d , where d is the embedding dimension. For 


example, if the elements are uniformly scattered on a smooth 
surface (d = 2), the probability of finding a neighbor within 
a given distance l scales as the area of a circle with this ra- 
dius, and hence C oc f 2 . If the elements are clustered, then 
there exists an interval of scales in which C(£) oc l “ <d . The 
power-law exponent a of the correlation integral is an empirical 
proxy for the correlation dimension. If a is constant across a 
substantial range of scales (one or more orders or magnitude), 
then the studied set demonstrates self-similar or self-affine den- 
sity variations involving voids and aggregates of various sizes. 
The discrepancy between the correlation dimension a and the 
embedding dimension d is a sensitive measure of this clustering 
behavior. 

The Cl technique can be generalized to the case when the 
studied set contains a mixture of elements of two distinct types 
(Kantz 1994). In our study, the two types of events are associated 
with positive and negative line-of-sight polarity magnetic events 
(origin or demise) embedded in space and in time. The events 
are described by a combination of spatial coordinates ri and 
timings f ; , i = 1, .., N , where /V = N+ + /V_ is the total number 
of positive and negative events. The subsets associated with OE 
and DE dynamics are treated independently. 

As the preparatory step, we compute the matrix L = { Lq } 
of pairwise distances Lq = r,- — r,j and the vector T = [Tq] 
of pairwise time intervals Tq = |f, — tj | between all the events, 
where i, j e [1, A]. We also define the binary vector P = { P, | , 
taking the value 1 (0) if the ith event is of a positive (negative) 
type. Next, we construct the inclusion operator 

r) = 0(1 - Lq)0( t - Tq ) (2) 

in which 0 is the Heaviside step function. Note that ^fq(i, t) 
equals 1 if the two events i an j lie within spatial distance l and 
their timing difference does not exceed r, and takes the value 
0 otherwise. Finally, we define the following set of generalized 
correlation integrals (GCIs): 

C+(£,T)= (A (3) 

c - (£ ’ T) = WY *) ( 4 ) 

C±(£, r) = 1 V Piil-PAVqil, r). (5) 

(A + )(A_) jy 1 1 

Each of the GCI functions (C+, C_, and C±) approximates 
the probability of finding a pair of events of a particular type 
(respectively, positive, negative, and bipolar) within the spa- 
tiotemporal neighborhood V = v,j r <x £ d t in which v,j is the 
volume of the (7-dimensional hypersphere (Figure 2). By con- 
struction, the autocorrelation integrals C+ and C_ characterize 
the internal structure of positive and negative subsets of events 
and are analogous to the usual Cl (Grassberger & Procaccia 
1983), whereas the cross-correlation integral C± describes in- 
teractions between positive and negative events in bipolar mag- 
netic pairs. The mathematical definition of C± is similar to the 
cross-correlation sum proposed by Kantz (1994) for evaluating 
an overlap of two chaotic attractors embedded in a reconstructed 
phase space of a nonlinear dynamical system. The primary 
difference is that here we deal with an essentially anisotropic 
metric space spanned by spatial and temporal coordinates of 
solar magnetic events. 


9 = 1 


otherwise. 
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Figure 2. Schematic illustration of the GCI technique for d = 2. Black (white) 
circles show positive (negative) events embedded in the d + 1 -dimensional space- 
time. The cylinder represents a spatiotemporal neighborhood centered around 
the event outlined by the square and defined by a particular choice of i and 
r. The outlined (positive) event has N+ = 7 unipolar and N± = 6 bipolar 
neighbors. The GCI calculation consists of an ensemble-averaging over all such 
neighborhoods centered successively at each of the events in the studied sets, 
for all relevant combinations of l and r . 

It is easy to show that randomly distributed events with 
statistically independent positions and occurrence times yield 
C oc U, VC e {C+, C_, C±}. Any departure from this “null” 
scaling regime indicates spatial and/or temporal correlations 
between the events. By adopting the scaling ansatz developed 
in the theory of nonequilibrium critical phenomena (Barabasi & 
Stanley 1995), we assume that C(£) and C(r) dependences can 
be factorized after a proper renormalization resulting in a set of 
spatial and temporal power laws: 


C+ OC t a *T P* 

(6) 

C _ oc£“-r /i - 

(7) 

C± oc l a± x p± . 

(8) 


polarity sets are generated by randomly sampling the same un- 
derlying Cantor set. They have a statistically similar shape, and 
all three correlation integrals (C+, C_ and C±) collapse onto 
a single curve, with the log-log slope approaching dp as ex- 
pected. In panel (b), the compared sets are created by sampling 
two independent runs of the Cantor algorithm. These sets have 
the same autocorrelation properties as in Figure 3(a), but must 
be uncorrelated with each other. Indeed, the measured a+ and 
a dimensions are close to the theoretical value dp, whereas 
a± approaches the embedding dimension d — 2 confirming mu- 
tual statistical independence of the two sets. The GCI method 
leads to this correct conclusion despite the complex structure of 
the compared sets generating false visual similarities at various 
scales. 

For the sake of brevity, the experimental part of this paper 
focuses on spatial scaling properties of magnetic events as 
reflected by a exponents computed at different r scales. The 
delay time is varied but treated as a parameter rather than 
as an independent scaling variable, with the main emphasis 
placed on the multiscale geometry of unipolar and bipolar 
magnetic structures as well as the mechanisms of their cross- 
scale interaction and organization. The analysis of ft exponents 
and the associated temporal clustering is left for future research. 

To estimate a dimensions, we evaluated log-log slopes of 
the corresponding GCI functions using a least-squares linear 
regression model built for several selected ranges of spatial 
scales. In addition to this approach yielding a discrete set of a 
values, we used the continuous scalogram technique introduced 
by Uritsky et al. (20 11) which enables classification of clustering 
regimes across the entire studied range of scales. The GCI 
scalograms are defined as 


a(£, x) = 


9 log C{i, r) 
9 log i 


(9) 


The spatial autocorrelation dimensions a+ and a_ character- 
ize geometric clustering of, respectively, positive and negative 
events; the events of either type are clustered if the correspond- 
ing dimension is below d. Temporal autocorrelation dimensions 
ft+ and ft - have the same meaning in the time domain, with 
ft < 1 reflecting clustering effects along the time axis (since the 
embedding dimension of time is 1). 

As opposed to 04, a_, ft + , and ft _ dimensions describing 
unipolar correlations within the sets of positive and negative 
magnetic events, the cross-correlation dimensions a± and ft± 
are intended to quantify mutual coupling between the two 
types of events. If negative and positive events are statistically 
independent in space and in time, then C± scales as oc r d r 
similarly to the autocorrelation integrals C+ and C_, and so 
a± = d and ft± = 1. For correlated sets of bipolar events, 
the cross-correlation dimensions take values below or above 
the embedding dimension depending on the sign of the bipolar 
correlation (positive or negative, correspondingly). 

Asymptotically, the cross-correlation dimensions are inde- 
pendent of the intrinsic clustering structure of the compared 
sets (Kantz 1994). In other words, if each of the unipolar sets 
is clustered but the sets are unrelated to each other, then we 
expect that a± » d and ft± ~ 1 for sufficiently large N and 
negligible finite-size effects (so that l and r are much smaller 
than the largest relevant spatial and temporal scales). Figure 3 
illustrates this property of the GCI statistics on static sets of 
bipolar events constructed by using the stochastic Cantor set 
algorithm (Turcotte 1989) with the theoretical fractal dimen- 
sion dp = log(3)/log(2). In panel (a), positive and negative 


where a e {a + , a± j and C e {C+, C_,C±). As in the 
case with the discrete GCI dimensions, the scalograms (9) were 
computed independently for the OE and DE populations of 
magnetic events. 

The statistics reported in the next section are for a single 
detection threshold B± = 150 G. Most of our conclusions hold 
for a wide range of B± within I 4 standard deviations above 
the mean magnetic field magnitude, with the exception of the 
cross-correlation dimension of the active region exhibiting a 
more subtle B± dependence as discussed in Section 3.3. 

3. RESULTS AND DISCUSSION 

3.1. Overview of Scaling Properties of the Magnetic Events 

Our analysis has shown that the cooperative behavior of the 
same and mixed polarity magnetic events depends strongly on 
their spatial and temporal separation as well as on the local 
solar activity level. It has also been found that bipolar events 
of either type (origin or demise) exhibit considerably different 
spatial correlation patterns compared to the corresponding type 
of the unipolar events. The difference is statistically significant 
and has a tendency to increase with coronal activity. 

Figure 4 presents the combined results of the GCI analy- 
sis of OE and DE activity detected in the quiet and active 
photospheric regions. The left panels show the GCI functions 
(Equations (3)-(5)) evaluated for two fixed timescales, r = 
0 minutes (dotted lines) and r = 10 minutes (solid lines). It can 
be seen that the shape of the autocorrelation integrals charac- 
terizing unipolar magnetic clusters changes dramatically with 
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Figure 3. Identifying auto- and cross-correlations in simulated clustered sets using the GCI approach. The studied sets (blue and green dots) are generated by the 
same (a) and different (b) realizations of the underlying stochastic Cantor set with fractal dimension log(3)/ log(2) « 1.158. Bottom panels (c) and (d) show auto- and 
cross-correlation integrals (Equations (3)-(5)) computed for the pairs of sets in (a) and (b), correspondingly. The autocorrelation integrals of each set are shown with 
the matching color (blue or green). The cross-correlation integrals are in red. The code of the stochastic Cantor set is courtesy of G. Uritsky. 


the delay time. For r = 10 min, both C+ and C_ exhibit a 
distinct concave crossover (3 2 C/di 2 > 0) in the vicinity of 
£ = Iq ~ 15 Mm. As we discuss below, this spatial scale 
can be associated with the lower size limit of photospheric 
supergranules. 

The right panels in the same figure show the dependence of 
the GCI dimensions measured below (dashed lines) and above 
(solid lines) the crossover scale Iq on the timescale r. The 
systematic decrease of the autocorrelation dimensions a+ and 
«_ is likely a result of the /’-mode contamination due to the 
5 minute Doppler oscillations leaking into the Zeeman signal 
and/or the granulation noise described by roughly the same 
timescale (DeForest et al. 2007). Each of these mechanisms 
can modulate the magnetic flux, causing the event detection 
algorithm to loose track of weak magnetic features with a 
field magnitude close to B±, and to rediscover them after one 
disturbance cycle. In either case, the increase of r beyond the 
5 minute timescale leads to better statistical representation of 
long-living photospheric structures. It is worth noticing that 
the cross-correlation integrals C± (shown in red) are generally 
steeper than C+ and C_, and unlike the autocorrelation integrals 
reveal no systematic dependence on r. 

Table 1 summarizes our measurements of the GCI dimensions 
for the quiet Sun as well as for the active region operating in the 
quiescent and flaring states. Sections 3.2 and 3.3 below provide 
a detailed discussion of the scaling properties characterizing 
these three regimes. 

3.2. Quiet Sun 

As can be seen from Figures 4(a) and (b) and Table 1, unipo- 
lar dimensions a+ and a_ of the quiet photosphere are sys- 
tematically below the value 2, revealing clustering effects at the 
levels of meso- and supergranulation. At distances l < Iq, these 


dimensions are in the range 0 to 1 for both origin and demise 
events. At larger spatial scales, the unipolar dimensions take 
values between 1 and 2. The observed scaling crossover uncov- 
ers two distinct clustering regimes characterizing the geometry 
of the photospheric flux below and above the supergranular 
scale Iq. The small-scale dimensions are most likely to rep- 
resent the tendency of the magnetic concentration regions to 
align along the boundaries of the photospheric supergranules. 
The non-integer a values indicate the presence of gaps and 
clusters of various sizes analogous to fractal clustering in a one- 
dimensional Cantor dust model (Turcotte 1989). The large-scale 
clustering at l > Iq described by a_, a+ e [1,2] reflects a more 
isotropic spatial arrangement of solar magnetic flux. At these 
scales, the flux should be organized into quasi-two-dimensional 
supergranular structures also containing multiscale voids and 
clusters. 

It is of interest to relate the obtained statistics with commonly 
discussed theoretical scenarios of photospheric supergranula- 
tion. The supergranular flow is known to contain a dominant 
horizontal component (e.g., Leighton et al. 1962; Hathaway 
et al. 2000). Its leading role in the formation of the magnetic 
network outlining supergranular cell boundaries has been inten- 
sively discussed (Hagenaar et al. 1997; Schrijver et al. 1997; 
Srikanth et al. 2000). In deterministic explanatory models, the 
photospheric network appears as a result of forced rearrange- 
ment of small-scale surface magnetic fields by a large-scale 
horizontal flow (Martin 1988; Lisle et al. 2000; Wang et al. 
1996). In stochastic models, the supergranulation network is 
generated by kinematic motion of elementary flux tubes trans- 
ported passively by the supergranular flow superposed with ran- 
dom displacements produced by granulation (Simon et al. 1995, 
2001), by mutual advection of thermal downflows taking place 
at random locations (Rast 2003), or as a self-organization of 
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Figure 4. Results of the GCI analysis of origin and demise activity in quiet and active photospheric regions. Left panels: C+, C . and C± plots for two timescales 
r = 0 (10) minutes shown with dotted (solid) lines. Three integer- valued log-log slopes are shown for reference, with a = 2 marking uncorrelated regimes. The 
boundaries of the two scaling ranges of interest associated with the spatial organization above and below the characteristic scale i o = 1 5 Mm are shown with vertical 
lines. Right panels: dependence of the dimensions a? + , a . and a± on r, in two ranges of spatial scales: l < lo (solid lines) and (. > lo (dashed lines). Color coding: 
blue (green) — unipolar C and a statistics of positive (negative) events, red — bipolar statistics of mixed polarity events. 


many interacting magnetic elements following a random walk 
pattern (Crouch et al. 2007). 

The deterministic models operate with a well-ordered super- 
granular flow which would result in integer a values both at 
small and large spatial scales, contrary to our observations. The 
second (stochastic) group of models seems to be in a qualita- 
tive agreement with the correlation analysis results. Our mea- 
surements show that the network is significantly fragmented 
at the supergranulation scale and exhibits an approximate co- 
alignment of magnetic flux elements with the lines of merging 


of the adjacent supergranular flow cells. The fractional values 
of a_ and o! + exponents at both small and large scales can be 
a footprint of a stochastic generating process, possibly related 
to the aggregation, fragmentation, and cancellation of randomly 
shuffled magnetic concentration regions as proposed by Crouch 
et al. (2007). 

To some extent, our estimates can be affected by the Doppler 
noise due to the flow motion at the supergranulation scale. 
In contrast to the short-scale noise sources mentioned in 
Section 3.1, the supergranulation flows are capable of 
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Table 1 

Estimated Values of GCI Dimensions of Origin (OE) and Demise (DE) Magnetic Events in the Quiet (Q) and Active (A) Regions 


Region 


Z <Z o 



l > Zo 



ct+ 

a- 

oe± 

a+ 

a- 

oe± 

Q 

1.66 it 0.20 

1.36 ±0.06 

OE, r = 0 minute 
2.65 d= 0.41 

1.64 ±0.01 

1.19 ±0.02 

1.86 ± 0.01 

A 

2.03 ± 0.27 

2.01 ±0.13 

2.06 ± 0.18 

1.46 ±0.01 

1.56 ±0.01 

2.24 ± 0.02 

Q 

1.64 ±0.13 

2.01 ±0.18 

DE, r = 0 minute 
3.64 ± 0.39 

1.64 ±0.01 

1.24 ±0.02 

2.47 ± 0.01 

A 

1.30 ±0.07 

1.77 ±0.05 

1.25 ± 0.06 

1.33 ± 0.01 

1.51 ±0.01 

1.95 ± 0.02 

Q 

0.43 ± 0.02 

0.72 ± 0.02 

OE, r = 10 minutes 
3.28 d= 0.13 

1.40 ±0.01 

1.01 ±0.01 

2.25 ± 0.01 

A 

0.78 ± 0.02 

0.55 ± 0.02 

2.02 ± 0.06 

1.23 ± 0.01 

1.43 ±0.00 

2.02 ± 0.01 

Q 

0.43 ± 0.01 

0.63 ± 0.02 

DE, r = 10 minutes 
3.53 d= 0.15 

1.46 ±0.01 

0.97 ± 0.00 

2.32 ± 0.01 

A 

0.86 ± 0.01 

0.68 ± 0.02 

1.84 ± 0.07 

1.21 ± 0.01 

1.41 ±0.00 

2.03 ± 0.01 

Q 

0.35 ± 0.02 

0.59 ± 0.02 

OE, r = 20 minutes 
3.10 ± 0.11 

1.38 ± 0.01 

0.96 ± 0.01 

2.26 ± 0.01 

A 

0.72 ± 0.01 

0.54 ± 0.02 

1.94 ± 0.06 

1.21 ± 0.01 

1.40 ±0.01 

2.00 ± 0.00 

Q 

0.39 ±0.01 

0.57 ± 0.02 

DE, r = 20 minutes 
3.31 d= 0.12 

1.38 ± 0.01 

0.95 ± 0.01 

2.33 ± 0.01 

A 

0.72 ± 0.01 

0.56 ± 0.02 

1.70 ± 0.06 

1.22 ±0.01 

1.37 ±0.00 

2.04 ± 0.01 


injecting additive signals into MDI magnetograms over a wide 
range of scales (see, e.g., Meunier et al. 2010 and references 
therein). However, we expect this mechanism to have a limited 
effect on our results. The rms values of the erroneous magnetic 
field caused by Doppler noise are known to increase with mean 
magnetic field. The error should be quite small for the magnetic 
field magnitudes comparable with the chosen detection thresh- 
old B ltl . The simulation study by Settele et al. (2002) predicts 
SOHO MDI rms errors to be less than 10 G for a relatively strong 
mean signal of 800 G subject to 5 minute velocity oscillations 
with amplitude up to ~300 ms -1 . Extrapolating their estimates 
to fi t h = 150 G yields an uncertainty of ~1 G which can be 
safely neglected. This conclusion applies to both quiet Sun and 
AR 10365 active region since we applied the same detection 
threshold on both sides of the solar equator. 

Furthermore, the supergranulation Doppler noise exhibits 
strong dependence on the heliocentric angle from the disk 
center. At large heliocentric angles, the line-of-sight motion 
is dominated by the horizontal component (100-300 m s -1 ; see, 
e.g., estimates by Welsch et al. 2012 based on local correlation 
tracking) which can lead to a sizable Doppler effect. At smaller 
angles, this motion is represented by the much smaller radial 
velocity component ranging from 7 to 30 m s' 1 (Williams & 
Pesnell 2011; Hathaway et al. 2002). Our tests have shown 
that all main features of the GCI statistics discussed in this 
section, including the crossover scale 1$, are clearly present in 
this central solar region where the supergranulation Doppler 
noise plays only a minor role. 

Unlike the autocorrelation (unipolar) dimensions, the a± 
dimensions of the quiet photosphere are systematically greater 
than 2, pointing at a consistently negative cross-correlation 
between opposite polarity events. To visualize the clustering 
geometries underlying the auto- and cross-correlation results, 
we constructed the unipolar and bipolar clustering maps defined, 
respectively, as M up = {r,- — r*. | P t = P ).} and Mb p = 
{ r, — r/ f | P; ^ P/.}. In essence, the maps are obtained by adjusting 
the origin of the coordinate system consequently with each of the 
N events, and superposing the results in a single scatter plot while 
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Figure 5. Unipolar and bipolar clustering maps of the quiet photospheric region, 
M up and Mb p , visualizing the difference in the geometry of the OE dynamics of 
the same (top) and opposite (bottom) magnetic polarity. 

(A color version of this figure is available in the online journal.) 

keeping track of the polarity of the events. As expected, the two 
maps are significantly different in shape (Figure 5). The dense 
core of the unipolar map represents the tendency of the magnetic 
events of the same sign to cluster into compact groups described 
by a± < 2. The empty central region of the bipolar map reflects 
negative correlations inside the bipolar pairs leading to a± > 2, 
with the probability of finding a neighbor of an opposite polarity 
increasing with the distance up to l ~ 30-40 Mm. 


unipolar 


bipolar 
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The observed anti-correlation of the bipolar events in the 
quiet-Sun network implies a significant imbalance of the ap- 
parent evolution of the magnetic flux in this region. Since the 
actual flux dynamics must be balanced (V ■ B = 0), our statistics 
manifest a geometric asymmetry of newly injected or cancel- 
ing bipoles precluding the detection of a matching pole. It is 
known that an asymmetric coalescence of a preexisting unre- 
solved unipolar flux into magnetic concentrations that are suf- 
ficiently strong to be detected can produce unbalanced OEs: if 
two or more like-polarity ends of previously unresolved flux 
tubes come together, the average field strength can overcome 
the detection threshold for that polarity only (Lamb et al. 2008). 
Likewise, asymmetric fragmentation can lead to unbalanced 
demise events. Such sign-dependent effects are consistent with 
our analysis of the quiet solar region. We find that the posi- 
tive events occurring in this region are by a factor of 1. 5-3.0 
(depending on r and Zf t h) more frequent than negative events, 
indicating that the positive magnetic flux was significantly more 
fragmented than the negative flux, which could naturally lead to 
the described artifacts. This interpretation is also supported by 
the fact that the dimension a + is systematically lower than a_ 
(see Table 1), suggesting that spatial arrangement of positive 
polarity events along supergranule boundaries is more sparse 
compared to that of negative events. 

In addition to being a potential source of false event detection, 
the sign-dependent coalescence-fragmentation dynamics imply 
that, on average, each negative event is magnetically connected 
with more than one positive event. The higher-order correlations 
associated with such magnetic structures are beyond the scope of 
applicability of the proposed second-order statistical measures 
focusing on pairwise bipolar interactions such as those in 
magnetic loops. Addressing such correlations in future studies 
would be a challenging but rather important task. 

It is also noteworthy that in spite of the geometric asymmetry 
of positive and negative polarity fluxes, the origin and demise 
GCI dimensions of either magnetic polarity obtained for the 
quiet Sun are roughly the same. This can be a manifestation 
of a steady-state balance between the rate of injection of the 
magnetic flux from the convection zone, and the rate of its 
dissipation via coronal heating. It is generally accepted that the 
network magnetic fields play an important role in heating the 
solar corona. The total X-ray brightness of an active region is 
closely connected with the overall length of neutral lines and 
the presence of strong magnetic shears in the neighboring area 
(Falconer et al. 1997). Active core fields of the network of the 
quiet corona can be the main drivers of bulk coronal heating and 
dissipation as evidenced by the correlation between the location 
of bright X-ray emission points and bipolar magnetic structures 
(Falconer et al. 1998). Placed in this context, the nearly identical 
origin and demise GCI statistics strongly suggest that the coronal 
heating in the quiet Sun can be a magnetically driven process 
in which the sources and the sinks of the magnetic free energy 
share the same long-range correlation structure. 

The possibility of an evolving but statistically stationary 
quiet-Sun network such as the one revealed by our analysis has 
been predicted earlier by Simon et al. (2001) based on kinematic 
modeling. Their results show that the total unsigned magnetic 
flux on the solar surface would decay within a few days due to 
the annihilation of oppositely directed fields brought together 
by granular and supergranular flows, unless it is continually 
replenished at an appropriate constant rate. This requirement 
was presumably fulfilled in the studied quiet solar region as 


suggested by the statistically identical scaling behavior of the 
OE and DE dynamics. 

3.3. Active Region NOAA 10365 

The correlation integral analysis of the active region leads to a 
rather different picture in which the origin and demise magnetic 
flux events follow distinct scaling laws, with a± » 2 for the 
former and a± < 2 for the latter events, see Figures 4(c) and (d). 

The a± & 2 condition describing the OE population suggests 
the absence of sizable spatial correlations in newly detected 
bipolar pairs. This can be a signature of turbulent dynamics 
in the convection zone or some other process leading to a 
highly fragmented polarity separation line. Such lines are often 
associated with a sustained shearing motion of the field leading 
to a complex mixture of opposite magnetic polarities and 
resulting in the nonpotentiality of the active region coronal field 
(Schrijver et al. 2005). Typically, the heating and flaring of 
coronal loops above an active region requires magnetic polarity 
inversion near at least one of the loop footpoints, and is greatly 
enhanced in the presence of a strong magnetic shear in the core 
magnetic field (Falconer et al. 1997). The disordered topology 
of the shear suggested by our GCI measurements may play a 
vital role in initiating magnetic reconnection above solar active 
regions (Antiochos 1998; Podgorny et al. 2007; Torok et al. 
2011 ). 

In agreement with the proposed interpretation, strong shear- 
ing flows were observed near the main polarity inversion line 
of NOAA 10365 (Chae et al. 2004). The region showed a sys- 
tematic buildup of magnetic complexity en route to a flare- 
productive unstable configuration. The onset of coronal activity 
of NOAA 10365 was preceded by a pronounced distortion of 
the shape of the main bipole separation line (Verbeeck et al. 
2011), an injection of the helicity in the polarity inversion re- 
gion (Chae et al. 2004), and a sharp increase of the Ising energy 
(Ahmed et al. 2010) of the bipolar magnetic flux (Verbeeck et al. 
2011). A multi-wavelength analysis confirms the formation of 
a highly twisted magnetic flux tube structure prior to M-class 
flares in the NOAA 10365 region. Up to 50% of the helicity in- 
jection rate was contributed by short-lived magnetic structures 
with characteristic times of 10-60 minutes (Tian et al. 2011). 
The precursive process outlined above could easily generate a 
rather complex photospheric field pattern with highly entan- 
gled magnetic polarities. In fact, the value a± rs 2 implies 
that the complexity of the magnetic network of the active 
region was close to the theoretical limit in which the emer- 
gence locations of oppositely directed flux elements are totally 
uncorrelated. 

In contrast, the demise events in NOAA 10365 are described 
by substantially lower a± values signaling positive correlations 
in bipolar magnetic pairs, possibly due to the cancellation of 
flux in reconnecting magnetic loops structures. This possibility 
is supported by the dependence of the cross-correlation dimen- 
sion on the level of coronal flaring activity as shown in Figure 6. 
The figure presents continuous scalograms a±(£, r) defined by 
Equation (9) characterizing the magnetic flux demise dynamics 
in flaring and non-flaring regimes of NOAA 10365 (panels (b) 
and (c), correspondingly) as compared to the demise dynam- 
ics in the quiet photosphere (panel (a)). Each scalogram is ac- 
companied by the r-dependence of small- and large-scale GCI 
dimensions plotted in the same format as the right column in 
Figure 4. 
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Figure 6. Comparison of the correlation structure of the demise events in the quiet photosphere (a), in the NOAA 10365 active region without significant flaring 
activity (b), and in the same active region producing M-class flares (c). Left column: GCI exponents vs. timescale r plotted in the same format as in Figure 4. Right 
column: continuous scalograms showing the dependence of the cross-correlation dimension a± on spatial and temporal separation between bipolar demise events. 
The red color prevailing in the flaring active region scalogram may be a signature of the magnetic flux cancellation. 


The level of flaring activity is clearly encoded in the shape of 
the scalograms shown in Figure 6. The flaring state of the active 
region is characterized by the condition a± < 2 (scalogram 
regions colored in red and yellow) for almost the entire studied 
range of l and r , confirming the presence of positively correlated 
bipolar demise dynamics at various spatial and temporal scales. 
The sign of the correlation changes from positive to slightly 
negative only in two narrow bands which could be associated 
with finite size effects imposed by the geometry of the active 
region. Even for very small delay times r ~ 1 min, there exist 
a distinct domain of l scales (~ 10-20 Mm) where the cross- 
correlation is positive. As was already mentioned in Section 3.2, 
the quiet Sun exhibits no stable positive bipolar correlations. 
The narrow small-scale i band with a± < 2 seen in the 
top right panel of Figure 6 can be an artifact reflecting an 
insufficient number of closely located bipolar pairs in the quiet- 
Sun network. The quiescent active region takes the intermediate 
position between the quiet Sun and the flaring active region, 
showing small but statistically significant positive correlations 


in selected i ranges. As in the case with the flaring NOAA 10365, 
the correlations describing nearly simultaneous demise events 
(r ss 0) in the non-flaring state of this active region are only 
detectable at short spatial distances. 

The observed tendencies suggest that the GCI dimensions 
of the active region indicative of positively correlated bipolar 
demise events are connected with the flaring activity. In other 
words, there is a coupling between the symmetric reduction 
of the photospheric flux and coronal dissipation. The corre- 
lated demise events could be produced by multiple bursts of 
magnetic reconnection in NOAA 10365 leading to spatially lo- 
calized symmetric reduction of oppositely signed line-of-sight 
magnetic fluxes. Indeed, the geometry of the chromospheric 
ribbons observed in the H a band unambiguously shows that 
the flaring activity of NOAA 10365 was triggered by magnetic 
reconnection at coronal heights below the twisted flux tube of 
positive helicity (Chandra et al. 2009). This is consistent with 
the results of 3D MHD simulations of the studied active region 
showing the formation of multiple current sheets in the vicinity 
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Figure 7. Examples of spatial arrangement of magnetic events in the quiet 
Sun (a) and in the active solar region in quiescent (b) and flaring (c) states 
consistent with our analysis results. Straight line segments indicate positively 
correlated demise events of opposite magnetic polarity presumably associated 
with magnetic flux cancellation events in the active region which could be 
involved in sympathetic solar eruptions (Torok et al. 2011). 

(A color version of this figure is available in the online journal.) 


of coronal X-lines above NOAA 10365 (Podgorny et al. 2007; 
Podgorny & Podgorny 2008). Each of the current sheets may 
well be responsible for an elementary flare. The reconnection- 
driven annihilation of the coronal field above NOAA 10365 
could lead to a detectable reduction of the photospheric flux in 
magnetically connected bipolar magnetic pairs, and contribute 
to the observed correlated occurrence of the demise events char- 
acterized by the condition a± < 2. To some extent, the bipolar 
GCI statistics can also represent simultaneous reduction of the 
magnetic flux in canceling bipolar pairs involved in bulk coronal 
heating around the active region (Schrijver & Title 2002). 

The r-dependence of a± dimensions of the active region is 
different than that of the quiet Sun, with positive correlations 
in disappearing bipolar pairs being the strongest at smallest 
r scales. This suggests that the characteristic timescale of the 
apparent flux cancellation can be comparable with the 1 minute 
image sampling time. The correlations also increase with the 
decrease of the spatial scale i and tend to be less pronounced 
at i > £(). If the association of these events with canceling 
magnetic flux elements is correct, then our observations indicate 
the dominant contribution of small- to medium size coronal 
loops in the reconnection process. 

Positively correlated unipolar demise events are reliably iden- 
tified for B t h exceeding ~2 standard deviations of the active 
region field but are not seen at lower detection thresholds. 
Therefore, they are the property of an intense magnetic field. For 
higher thresholds, the cross-correlations monotonously increase 
with B t h suggesting the involvement of non-Gaussian magnetic 
field fluctuations, such as those associated with intermittent tur- 
bulence below the photosphere (Abramenko 2005) or transient 
multiscale dissipation in the corona (Uritsky et al. 2007). This 
explanation is in line with earlier observations showing that the 
preferred locations of bright coronal points correlate with re- 
gions of intense unsigned magnetic flux (see Benevolenskaya 
2010 and references therein), although the physical mechanism 
of this dependence, likely driven by topological transitions in 
a complex multiscale magnetic network, remains to be under- 
stood. 


4. CONCLUSIONS 

Our findings demonstrate that the developed data analysis 
approach based on the GCI algorithms and direct spatiotem- 
poral tracking of magnetic elements is capable of addressing 
various aspects of the dynamic of solar photospheric network, 
including the geometry of supergranular structures, emergence 
of the nonpotential flux, and cancellation dynamics of bipolar 
magnetic pairs. 


The results obtained reveal significantly different multiscale 
dynamics of the photospheric magnetic flux in quiet and active 
solar regions as represented by the studied sample of SOHO 
MDI images. The quiet Sun exhibits robust unipolar clustering 
of both OE and DE events over a broad range of photospheric 
distances corresponding to meso- and supergranulation, with a 
small-scale quasi one-dimensional clustering along supergran- 
ular boundaries and a more isotropic clustering at the level of 
supergranules. These clustering regimes place important con- 
straints on the topology of the photospheric network and spa- 
tially distributed storage of free magnetic energy in the quiet Sun 
as illustrated by Figure 7(a). No detectable cross-correlations 
in bipolar magnetic pairs which could be related with unsta- 
ble magnetic loops were found in the quiet photospheric re- 
gion. NOAA 10365 solar active region exhibits signatures of 
non-potential magnetic structures with a topologically complex 
polarity separation line. 

Flaring activity in the active region is accompanied by posi- 
tively correlated magnetic demise events which can be caused by 
magnetic reconnection. A possible underlying geometry which 
would lead to the observed GCI signatures of the active region is 
a sheared loop arcade with an irregular polarity inversion line as 
depicted in Figures 7(b) and (c). The topological complexity of 
this structure can be instrumental because simple shearing of an 
arcade is typically insufficient for any significant manifestation 
of solar activity such as flares or prominences (Dahlburg et al. 
1991). The statistical results obtained for AR 103675 speak in 
favor of so-called sympathetic eruptions in this complex active 
region. Such eruptions may occur within a relatively short pe- 
riod of time in spatially separated flux ropes due to a variety of 
MHD-scale triggers destabilizing the coronal field (Torok et al. 
2011 ). 

The proposed methodology can be helpful in resolving a 
broader set of theoretical and practical problems associated with 
solar magnetic complexity and its impact on solar activity. For 
example, is there a relationship between statistical properties 
of fractality /complexity in the photosphere and corona? Some 
studies provide a definitive answer (e.g., Conlon et al. 2010; 
Falconer et al. 2009; Abramenko & Yurchyshyn 2010a, 2010b), 
but some others argue against such a possibility (Dimitropoulou 
et al. 2009). If the complexity linking the two solar structures 
exists, then one can expect the GCI signatures of the photosphere 
and corona to vary in a consistent way. 

It is also not clear how the fundamental plasma processes 
responsible for the formation and evolution of the magnetic 
network organize themselves across such a broad range of 
scales. Accurate evaluation of GCI dimensions for statistically 
representative sets of solar regions and conditions may provide 
important clues on the physics of the underlying process; 
for instance, by validating certain types of turbulent cascades 
(Muller & Biskamp 2000). The proposed methodology can 
also help to better understand how the scaling properties of 
the quiet and active Sun networks relate to each other, and 
how these properties reflect the formation of unstable magnetic 
topologies leading to plasma heating and dissipation. The 
relationship between the temporal and spatial scaling laws 
governing multiscale geometry and dynamics of the solar 
magnetic field is another important subject for future research. 
Scaling characteristics of sheared magnetic structures generated 
by oppositely directed line-of-sight flux elements, and the 
contribution of these structures to the coronal dissipation in quiet 
and active solar regions are also largely unexplored and require 
careful consideration. To clarify the underlying physics of the 
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observed effects, it would be essential to combine our statistical 
tools with conjugate observations of coronal activity as well as 
solar tomography (Davila 1994) and other three-dimensional 
reconstruction techniques. 

The developed approach can also be used for improving the 
forecasting of flaring activity through identification of topolog- 
ically complex unstable magnetic configurations. Spatial mag- 
netic complexity of an active region as measured by the power- 
law decay of the Fourier spectrum of the longitudinal field was 
shown to correlate with the flare productivity of the region 
(Abramenko 2005). This dependence suggests that active re- 
gions with more intermittent magnetic structure tend to produce 
higher energy flaring events. Our methodology offers a more 
comprehensive way of measuring the magnetic intermittency 
characterizing each magnetic polarity (based on a+ and a_), as 
well as bipolar magnetic structures ( a ± ). Analysis of these cor- 
relation measures, as well as the temporal GCI dimensions not 
considered here, may help us to better understand photospheric 
drivers of bulk coronal heating and local topological precursors 
of bright coronal points and flaring activity across the entire 
range of the involved spatial and temporal scales. 
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